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ABSTRACT 

Using multi-epoch observations of the Stripe 82 region done by Sloan Digital Sky Survey, we measure 
precise statistical proper motions of the stars in the Sagittarius stellar stream. The multi-band pho- 
tometry and SDSS radial velocities allow us to efficiently select Sgr members and thus enhance the 
proper motion precision to ~ 0.1 masyr^ 1 . We measure separately the proper motion of a photomet- 
rically selected sample of the main sequence turn-off stars, as well as of a spectroscopically selected 
Sgr giants. The d ata allow us to determ ine the proper motion separately for the two Sgr streams in 
the South found in lKoposov et al.l ()2012h . Together with the precise velocities from SDSS, our proper 
motion provide exquisite constraints of the 3-D motions of the stars in the Sgr streams. 
Keywords: Galaxy: halo, stars: kinematics, methods: statistical, surveys 



1. INTRODUCTION 

The disintegrating Sagittarius (Sgr) dwarf galaxy re- 
mains a riddle, wrapped in a mystery, inside an enigma. 

Large scale photometric surveys, such as the Two Mi- 
cron All-Sky Survey (2MASS) and the Sloan Digital Sky 
Survey (SDSS) have now revealed the structure of the 
tidal tails of the Sgr over more than 2ir radians o n the 
Sky (jMaiewski et al.|[200l iBelokurov et al.l [20061 ) . By 
tallying all the stellar debris in the streams and remnant, 
we now know that the progenitor galaxy had a luminos- 
ity of ~ 1O 8 L0, c omparable to the present day Sma ll 
Magellanic Cloud (|Niederste-Ostholt et al]|2010l I2012D . 
The ingestation of such a large progenitor, together with 
its dismantling under the actions of the Galactic tides, 
can provide us with a wealth of information about both 
the Galaxy and the Sgr, if we can only decode it. 

Radial velocities, and sometimes metallicitics and 
chemical abundances, are now known for many hun- 
dreds of stars in the Sgr tails (e.g.. iMaiewski et al.l 
[200llMonaco et al.ll2007tlChou et al.ll2007ll2010D . There 
are al so ~ 10 globular clusters associated with the Sgr 
tails (|Law fc Maiewskil l2010b[ ). This rich mosaic of 
positions and velocities of Sgr tracers has proved sur- 
prisingly difficult to understand. Although there is no 
shor tage of Sgr disruption models in the literature (see 
e.g.. lHelmi|[200l lLaw et al.|[200l I Johnston et~all [20051: 
iFellhauer et al.l I2006D . thev all have significant short- 
comings, and fail to reproduce a substantial portion of 
the datasets. The most su ccessful recent attempt is 
by lLaw fc Majewskil (|2010a[ ). though they do not ex- 
plain the st riking two stream morphology seen in the 
SPSS data (IBelokurov et al.|[200l IKoposov etail [20121: 
ISlater et al.ll2013l ). Additionally they advocate the use of 
a triaxial halo for the Galaxy with minor axis contained 
with the Gal actic plane, which is unatt ractive on other 
grounds (e.g. [Kuijkcn fc Tremainelll994T ). Given the im- 
passe, it is natural to look to proper motions of the Sgr 
stream as so as to obtain a clearer picture of its space 
motion. 
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iCarlin et al.l <|2012f ) have provided the first measure- 
ments of the proper motion of the Sgr trailing tail. 
They took advantage of archival photographic plate data 
in some of Kapteyn's Selected Areas which provides a 
90 year baseline. They derive proper motions for four 
40' x 40' fields covering locations on the trailing tail be- 
tween 70° and 130° from the Sgr core. However, the 
number of stars in each field remains modest (~ 15 — 55), 
and so the precision of the proper motion measurement 
is still quite low (~ 0.2 — 0.7 mas yr -1 ). 

Here, we will pursue a different tack to obtain proper 
motions of the trailing stream in roughly the same area 
of sky. As part of a project to detect supernovae, the 
Sloan Digital Sky Survey scanned a ~ 290 square de- 
gree regi on on the Celestial Eq uator, known as Stripe 
82 (e.g., lAbazaiian et al.l 12009ft . Proper motions 
can be derived by matching objects between the ~ 80 
epochs (jBramich et al.ll2008h obtained over time period 
of ~ 7 years, whilst the co-added optical data is roughly 
2 magnitudes deeper than a single epoch SDSS measure- 
ment. Although the baseline is small so the precision of a 
measurement of proper motion of a single star is still low, 
we can take advantage of the large number of Sgr tracers 
to get a high precision (~ 0.1 masyr -1 ) measurement for 
the proper motion of the ensemble. 

The paper is arranged as follows. Section 2 describes 
the extraction of proper motions for stars from the Stripe 
82 data using background quasars to provide an absolute 
reference frame. Section 3 shows how to identify the Sgr 
stars in Stripe 82, where they occupy a distinctive niche 
in magnitude and radial velocity space. Section 4 dis- 
cusses our modelling of the proper motions using both 
photometric and spectroscopic samples. We extract the 
proper moti on for both the bright and faint Sgr streams 
identified bv lKoposov et al.l(|2012fl . In section 5, we com- 
pare our proper m otions both with the earlier work of 
ICarlin et al.l (|2012[ ) and with the simulation data. 

2. STRIPE 82 PROPER MOTION DETERMINATION 

Stripe 82 has already been the subject of numerous 
studies. Its multi-epoch and multi-band imaging al- 
lows study of the variable sky and iden tification of many 
kinds of transient phenomena (see e.g..lSesar et al.ll2007l: 
iBecker et a l. 2008; Kowals kTeTall 120091: iWatkins et"atl 
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2009). The Stripe 82 dataset has also been used to de- 
rive proper motions bv lBramich et all (|2008| ). This light- 
motion catalogue was subsequen tly exploited to bui ld re- 
duced proper motion diagrams (jVidrih et alJl2"o"07T ) and 
analyse kine matical properties o f Galact ic disk and halo 
populations ([Smith et al.ll2009al lb1. [2012] ). 

Even so, the proper m otion measurements pioneered 
bv lBramich et al.l ([20081) can be improved. For the bulk 
proper motion of Sgr, we are interested in the statistical 
properties of a large ensemble of faint tracers, so proper 
motions with the smallest possible systematic errors are 
highly desirable. The original catalogue by Bramich does 
not provide proper motions for stars fainter than r ~ 
20.5 and is known to have some noticeable systematics. 
Since the Sgr stream has a very la rge number of tracers 
in Stripe 82 ([Watkins et al.l 120091 ) . we do not require a 
proper motion measurement for every star, but rather 
need small systematic errors and well understood error- 
bars for an ensemble. For this purpose, it makes sense to 
measure the proper motions relative to quasars (QSOs). 

Stripe 82 has a number of both spectroscopically 
and photometrically identified QSOs. In this work, we 
have used the c atalog ue of spectroscopic QSOs from 
ISchneider et al.l ([20101) an d the sample o f phot ometri- 
cally identified QSOs from iRichards et al.l ([2009D to ex- 
tract denizens of Stripe 82. The purity of the spectro- 
scopic catalogue of QSOs is guaranteed - all the objects 
are QSOs and must have zero proper motion. How- 
ever, the photometric catalogue is known to have some 
contamination by stars. In order to minimize contam- 
inants, we use the c ut good> 1, as recommended by 
IRichards et al.l ([2009! ) . This guarantees a small stellar 
contamination, certainly <5%. 



2.1. Relative Proper Motions 

Given a sample of QSOs each with zero proper motion, 
then for each star in the vicinity of the QSO, we may 
determine proper motion relative to the quasar. As an 
input ca talog for the stars , we took the Stripe 82 co-add 
dataset (jAnnis et al.l20li[) , from which we select primary 
objects, classified by the SDSS pipeline as stars. The 
individual source detections are taken from the Stripe 
82 po rtion of the SDSS DR7 database (|Q'Mullane et al.l 
120051 ) using only those fields having acceptable and 
good data quality flags. Matching co-added sources to 
detections at individual epochs is done with the 0.5 arc- 
sec radius using the Q3C module fo r the PostgreSQL 
database ([Koposov fc Bartvm ov 2006). This procedure 
makes the catalogue incomplete for high proper motion 
objects (with proper motions > 100 mas yr _1 ), but we 
are not interested in such objects in our current study. 

Then, for each pair (star, QSO) observed multiple 
times by SDSS within one field, we analyse the posi- 
tional offsets and errors. For the right ascension, these 
are defined via 
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with similar equations for the declination. The model for 



where /i is the proper motion of the star, t is the date 
of the observation, / is the fraction of outliers, and a is 
the scatter around the linear relation, whilst R(A) is the 
rectangular function to account for the outliers. 

The resulting likelihood is then minimized with respect 
to the 4 parameters /i, Ao, /, o with the error-bars de- 
termined from the Hessian at the minimum. Repeating 
this procedure for the offsets in declination gives us the 
proper motions and their errors, [i a , ct^q, fis, v^^g for all 
the sources with a spectroscopic or a photometric QSO 
nearby. 

2.2. Systematic Errors in the Proper Motions 

After performing the computation of the proper mo- 
tion for individual stars, and before trying to measure 
the statistical proper motions for ensembles of stars, it is 
important to check for the presence of possible systemat- 
ics, as well as to examine the accuracy of the error bars. 
Fig. Q] presents such an assessment. In this paper, we are 
focusing on the particular part of Stripe 82, which in- 
tersects with the Sgr stream. As systematic effects may 
depend on the right ascension, Fig.[T]uses only the proper 
motions in the right ascension range 20° < a < 50°. 

The left panel of Fig. Q] shows the histogram of proper 
motions of the spectroscopic QSOs measured relative to 
the photometric QSOs and normalized by the error bar 
provided by our fitting procedure. The red overplotted 
curve shows a Gaussian with the center at zero and unity 
dispersion. The excellent match of the histogram with 
the Gaussian curve shows that that the proper motions 
do not possess noticeable systematic offsets, and that the 
error bars on the proper motions are a faithful descrip- 
tion of the precision. It is also quite clear from the plot 
that this is not true f or proper motion measurements by 
IBramich et al.l (|2008[) , which have noticeable systematics 
and error underestimation. The middle panel of Fig. Q] 
shows the the proper motion of the spectroscopic QSOs 
relative to the photometric QSOs versus the g — r color 
difference. The points with error bars show the median 
proper motion in bins of g—r, while the error bars are 1.48 
times the median absolute deviation of proper motions 
within a corresponding g — r bin. Since all the points lie 
within 1 er of zero, we conclude that the color-dependent 
terms in the proper motions are negligible down to the 
precision 0.1—0.2 mas yr _1 . This holds true at least 
within the color range —0.2 < g — r < 1, which is the 
color-range applicable to 99% percent of the photomet- 
ric QSOs. Later we will see that the proper motion of 
the Sgr stream determined from the photometric sample 
with g—r ~ 0.3, g—i ~ 0.3 and the spectroscopic sample 
with g—r ~ 0.5, g—i ~ 0.75 agree each other, further con- 
firming the small level of color-related systematic errors. 
And last, the right panel of Fig. [T] shows the precision 
of our proper motion measurements as a function of r- 
band magnitude. There is a considerable scatter, but the 
median proper motion precisi on of ~ 2 mas yr -1 is sig- 
nificantly better than that of IBramich et al.l ([2008D for 
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Figure 1. Left panels: Normalized histogram of proper motions of spectroscopic QSOs measured relative to the photometric QSOs 
normalized by the error bar provided by the modeling. Red lines are G aussians with zero m ean and unity dispersion. Blue dashed curves 
show the histograms of the proper motions of spectroscopic QSOs from Bramich ct al. (2008). Central panels: The median proper motion 
of spectroscopic QSOs relative to the photometric QSOs versus their color- difference. Right panels: The measured error bars on the the 
proper motions as a function of the r band magnitu de. The red lin e s show the median values of our measurement errors. The blue dashed 
lines show the median of proper motion errors from Bramich ct al. (2008). 



Table 1 

Radial velocities and dispersions 
of Sgr stars along the stream as 
traced by SDSS. 



A 


Vgsr 


CTV 


deg 


krns" 1 


kms~ 1 


87.5 


-85.6 ± 2.1 


13.7 ±2.1 


92.5 


-100.5 ±1.4 


14.8 ± 1.5 


97.5 


-106.0 ±1.5 


11.1 ± 1.6 


102.5 


-121.0 ±1.0 


13.0 ± 1.0 


107.5 


-132.5 ±1.2 


13.8 ± 1.3 


112.5 


-140.5 ±1.6 


19.5 ± 1.4 


117.5 


-146.7 ± 1.5 


16.7 ± 1.4 


125.0 


-154.7 ±2.0 


19.9 ± 2.8 


135.0 


-156.8 ±2.8 


27.2 ±3.3 



large range of magnitudes. 

3. THE SAGITTARIUS STREAM 

The Sgr stream in the southern Galac tic hemisphere 
is kno wn to have a complex structure. IKoposov et al.l 
(2012) used main-sequence turn-off (MSTO) stars ex- 
tracted from SDSS Data Release 8 (DR8) to demonstrate 
the existence of two streams - a thicker, brighter stream 
and a thinner, fainter stream displaced by ~ 10°. The 
brighter stream has multiple turn-offs as well as a promi- 
nent red clum p , whe reas the fainter stream does not. 
IKoposov et al.l ()2012[ ) argued that the brighter stream 
was composed of more than one stellar population, in- 
cluding some metal-rich stars, whereas the fainter stream 
is dominated by a single metal-poor population. 

Here, we are primarily interested in the intersection 
of the Sgr stream with the SDSS Stripe 82. Fig. CU 
shows the density of MSTO stars extracted via the cuts 
0.25 < g-i < 0.35 and 19.8 < r < 22.5 in the south- 
ern hemisphere, together with Stripe 82 demarcated by 
the red lines. Again, two distinct streams structures are 
clearly visible - a brighter one crossing the Equator and 
Stripe 82 at a right ascension of ~ 35° and a dimmer one 
crossing the Equator at a right ascension of ~ 15°. The 
magnitude distribution of MSTO stars along Stripe 82 



is shown in the middle panel of Fig. O The Sgr streams 
are clearly evident, with the brighter stream visible at 
a ~ 35° and fainter one at a ~ 15°. Since the distances 
to the streams are not constant with a, and since Stripe 
82 crosses the stream at an angle, we observe a clear 
distance gradient with right ascension. 

The right panel of Fig. [5] shows the radial velocity as a 
function of right ascension for giant stars, extracted by 
the cut log(p) < 4 using the SDSS spectroscopic mea- 
surements. This time not only are stars within Stripe 82 
included, but also those satisfying |5| < 20° and lying 
ne ar the stream ( —15° < B < 5°). As already found 
bv lWatkins et al. (2009), the Sgr stream is visible at the 
Vgsr 150 km s" 1 B The fi gure also shows the vari- 
ation of radial velocities as the stream stars return from 
the apocenter of the Sgr orbit. As they will be needed 
later, we extract the radial velocities along the stream, 
together with the velocity dispersion, by fitting a Gaus- 
sian to the Sgr signal. Although Figure [5] shows the 
radial velocities as a function of right ascension, it is im- 
portant to realise that we perform these fits in bins along 
A, which is the angle along the stream from the remnant, 
measured positive in t he trailing direction, as defined in 
iMaiewski et al.l (|2003f) V The resulting measurements as 
a function of A are shown in Table [TJ 

Importantly, Fig. [5] demonstrates that the Sgr stream 
stars are visible as a distinct bright feature in magnitude 
and velocity space. Therefore, this information can be 
used to select Sgr member stars and statistically measure 
their properties such as proper motion. 

4. MODELING PROPER MOTIONS 

As demonstrated above, in the Southern Galactic 
hemisphere, it is possible to achieve clean separation of 
the Sgr trailing tail stars and the Galactic foreground in 

3 Here and th r oughout the pap e r, we use ^5^=23 5 km -1 
(IBovv et al.l [20091 ; IReid et al.l [2fjQ9l ; IDeason et all \20TH) and a 
solar peculiar velocity o f (U,V,W)=(-8.5,13.38,6.49) km s" 1 from 
HGoskunoelu et al.ll20TlT) 
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Figure 2. Left panel: Density of MSTO stars (0.25 < g-i < 0.35, 19.8 < r < 22) from SDSS DR8 in the southern Galactic hemisphere. The 
SDSS Stripe 82 region is delineated by the red line. Four horizontal error-bars show the range of right ascensions used to perform proper 
motion measurements in Stripe 82; red error-bars correspond to photometric samples, blue error-bars correspond to the spectroscopic sample. 
Center panel: The density of MSTO stars as a function of r band magnitude and right ascension illustrating a distance gradient along the 
stream. Right panel Radial velocities of stars near the Stripe 82. These are SDSS measurements of stars with log(g) < 4 and located near the 
Sgr orbit — 15° < B < 5° and near Stripe 82 \S\ < 20° . The Sgr stream is obvious at radial velocities of —200 km s _1 < Vgsr < —100 kms" 1 . 
The change of radial velocity along the stream is also quite prominent. 
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Figure 3. The bright Sgr stream component at 25° < a < 40° Left and Middle panels: Greyscale shows the 2D distribution of proper 
motions and magnitudes, while the red contours show the total error- deconvolved Gaussian mixture model. The component of the model 
corresponding to the stream is shown in blue. Right panel: ID projection onto the apparent magnitude axis of the left panel. The histogram 
of the data is shown in black, red curve shows the Gaussian mixture model of the luminosity function, whilst the model of the stream 
contribution is shown in blue. 



both apparent magnitude and the radial velocity space. 
However, even for high-confidence stream members, the 
uncertainty of individual proper motion measurements 
(~ 2 — 3 mas yr _1 ; see Section [2]) is comparable or higher 
than the expected tangential velocity of the stream (~ 
1 — 3 mas vr~ 1 : lLaw fc Majewskill2?IT0ah . It is therefore 
crucial to combine the signal from as man y stream mem- 
bers a s possible to beat down the noise. IKoposov et al.1 
(2010) have shown that, for the regions of apparent 
magnitude and radial velocity space dominated by the 
stream, an accurate measurement of the systemic proper 



motion of an ensemble of stars belonging to the stream 
can be obtained through simple background subtraction. 

Alternatively, the overall stellar distribution in the 
space of observables can be modeled, yielding the di- 
rect contributions of the Galaxy and the stream. If such 
models can be cast in the 3D space of proper motion and 
magnitude or radial velocity, it is feasible that a supe- 
rior measurement of the proper motion can be achieved 
as the contribution of the background to both system- 
atic and random noise will be reduced. Unfortunately, 
adequate analytical models of the Sgr stream and the 
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Figure 4. Left and Middle panels: Greyscale shows 2D histograms of proper motions and the GSR radial velocities corrected for the 
stream's radial velocity gradient for all stars with SDSS spectra and log(g) <4 in Stripe 82. Red contours show the error-deconvolved 
Gaussian mixture model of the data. Right panel: The histogram of the radial velocities. The Gaussian mixture model of the radial 
velocities is shown by the red curve and the stream component in blue. 



Table 2 

Proper motions measurements 



Field 


deg 


a 2 
deg 


/i a cos(S) h 
mas / yr 


mas/yr 


/A, a cos(5) 

mas/yr 


mas / yr 


Hi cos(b) b 
mas/yr 


mas /yr 


°7M 
mas / yr 


mas/yr 


FP1 


5.00 


22.00 


0.05 


-2.51 


0.11 


0.11 


0.10 


-2.51 


0.11 


0.12 


FP2 


25.00 


40.00 


-0.50 


-2.24 


0.09 


0.09 


0.81 


-2.15 


0.09 


0.11 


FP3 


42.00 


52.00 


0.19 


-2.36 


0.14 


0.14 


1.87 


-1.45 


0.11 


0.15 


FS4 a 


14.00 


50.00 


-0.42 


-2.65 


0.11 


0.13 


1.07 


-2.46 


0.11 


0.11 


Note. 


- ai 


CK2 columns denote the edges 


of the boxes 


in right 


ascension in 


Stripe 82 used to perform the 



proper motion measurements 



b not corrected for the Solar reflex motion 
a Spectroscopic sample 



Table 3 

Positions and velocities of the Sgr stream. 



Field 


A 


B 


X 


Y 


Z 


Distance 


U 


o"C7 


V 


cy 


W 


aw 




deg 


deg 


kpc 


kpc 


kpc 


kpc 


krns" 1 


krns -1 


km s — 1 


km s~ 


1 krns -1 


km s — 1 


FP1 


89.5 


10.0 


-15.0 


9.6 


-22.6 


25.4 


194 


14 


-46 


14 


23 


8 


FP2 


106.1 


0.6 


-23.5 


5.2 


-24.4 


29.1 


290 


12 


33 


14 


-20 


9 


FP3 


118.6 


-6.7 


-31.4 


0.4 


-25.0 


33.9 


270 


17 


-57 


24 


-48 


16 


FS4 a 


105.6 


0.8 


-23.2 


5.3 


-24.3 


28.9 


308 


15 


-16 


18 


-40 


11 



Note. — A, B, X, Y, Z correspond to the centers of the fields where the proper motions are measured. 
a Spectroscopic Sample 



Galaxy are not readily available. Therefore, we choose 
to approximate these distributions by a sum of multi- 
dimensional Gaussians. This so-called Gaussian mixture 
is a well known semi-parametric tec hnique widely used 
to mo del multi-dimensional datasets (McLach lan fc Peell 
2000). Gaussian mixtures have several key properties 
that make the model-fitting straightforward and fast. 
First, the uncertainties associated with the measure- 



ments (as well as the missing data) are naturally incor- 
porated into the model. Secondly, there exists a guaran- 
teed fast convergence proce dure - the Expectation Min- 
imisation (EM) algorithm (jDempster et al. 1 119771 ) . In 
this work, we have used the extreme-deconvolution pack- 
ag e, the open- s ource Gaussian mixture implementation 
bv lBovv et all (|20T1 . 
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Figure 5. Sgr streams in the South with the proper motion 
vectors ovcrplotted. The red vectors indicate the measurements 
performed using the photometrically selected MSTO stars at 
three different locations along the Stripe. Blue vector is for 
the spectroscopic sample. The measured error-bars of proper 
motions are shown in pink. The proper motion vectors have 
been corrected for th e solar motion, assuming the distances from 
IKoposo v et al\ H2012I) and Vlsr = 235 km s and peculiar veloc- 
ity from Coskunoglu ct al. (2011). 

4.1. Photometric sample 

As evident from the dissection of the SDSS dataset 
shown in Fig. [5J main sequence turn-off stars (MSTO) 
are the most numerous Sgr stream specimens available 
to carry out the proper motion analysis. The typical 
density of the MSTO stars in the stream is 150 stars per 
square degree, compared to the foreground density of 
around 50. Stripe 82 slices through both the bright and 
the faint streams at an angle, resulting in a slightly tilted 
bi-modal distribution in the plane of right ascension and 
r band magnitude. A faint Eastern wing to the bright 
Stream can also be discerned at higher right ascension 
(see middle panel of Fig. [2]). 

Motivated by the evolution of the distance and the 
structure of the stream along the Stripe, we split the 
MSTO sample into three parts: the faint stream at 
5° < a < 22°, the bright stream at 25° < a < 40°, and 
the Eastern wing of the bright stream at 42° < a < 52° 
(we label those fields as FP1, FP2, FP3 respectively; they 
are shown by red error-bars on Fig. [5]). In each of the 
three right ascension bins, the overall 3D distribution 
of MSTO stars in the space of (/i a ,/X5,r) was modeled 
with a mixture of 5 Gaussians, The initial guess for the 
free param eters was obtaine d by running the K-means 
algorithm (M acQueen I Il967t) . The choice of number 
of Gaussians (iV gau ) was mot ivated by looking at the 
cross- validated log-likelihood (jArlot fc Celissd 120101 ) as 
a function of -/Vg au . This initially rises as a function of 
A'gau, peaks at N gau ~ 5 and stays roughly constant for 
Ngau >5, demonstrating the goodness of fit and the ab- 
sence of overfitting. We have also made sure that the 
objects in the Sgr stream are represented by a single 



Gaussian. When applying the extreme deconvolution, 
we force the covariance matrix for the Sgr component to 
be diagonal, e.g. assuming zero covariance between the 
proper motion and the apparent magnitude. The proper 
motions for the stream stars were not corrected for the 
solar reflex motion. The errors in the proper motion 
measurement, i.e. the uncertainties in determining the 
centre of the Gaussian representing the Sgr Stream, were 
determined either from bootstrap procedure or from the 
Hessian matrix of the likelihood function. 

Fig. [3] shows the data analysed as well as the best 
fit Gaussian mixture model of the main Sgr stream 
25° < a < 40°. It is reassuring to see that the Gaussian 
mixture model was able to describe the data distribution 
adequately. The resulting measurements of the proper 
motion together with the uncertainties are given in Ta- 
ble [21 Among ~ 7000 stars analyzed, according to the 
model, ~ 4000 stars belong to Sgr, while ~ 3000 stars 
belong to the background/foreground population. 

4.2. Spectroscopic sample 

The number of members of the Sgr stream with SDSS 
spectroscopic measurements is significantly smaller than 
the number of MSTO stars used in the previous section. 
Despite this, knowledge of the radial velocities and sur- 
face gravities allows us to have much purer samples of the 
Sgr stream members. In this section, we perform a com- 
plementary measurement of the stream's proper motion 
using spectroscopic members only. 

To this end, we select all stars with spectra in Stripe 82 
lying at 14° < a < 50° ( we label this area FS4 and show 
it with the blue error-bar on Fig. [2]) and classified by the 
SDSS spectroscopic pipeline as giants log(g) < 4. De- 
spite very wide range of selected right ascensions, most 
of the Sgr members with spectroscopy in the Stripe 82 
region located at the center of the bright stream, at 
a ~ 30°. According to the right panel of Fig. [21 the 
Sgr stream's radial velocity changes along the Stripe. 
Therefore, to ease the modeling, we subtract the vari- 
ation of the radial velocity centroid from the data using 
the measurements presented in Table [TJ The resulting 
radial velocity is approximately constant as a function of 
RA, and the measurements of proper motions and radial 
velocities iu, a , fi s , V G sr = V GSR - Kiodd,GSR(A) can now 
be represented by a mixture of Gaussians. 

We run the extreme deconvolution on the sample of 
^1500 stars and find that 3 Gaussian components are 
sufficient to describe the dataset. We follow the strategy 
outlined in Section 14.11 with the difference that the co- 
variance matrix of the Gaussian representing the stream 
component is now set free. Fig. [4] shows the density dis- 
tribution of the data, together with the best-fit Gaus- 
sian mixture model for the Stripe 82 stars with spectra. 
By comparing the grey-scale density with the red con- 
tours, we can confirm that the model has captured the 
properties of the dataset reasonably well. The last line 
of Table [21 reports the values of the proper motion for 
the spectroscopic sample, and confirms that the mea- 
surements for the two independent stream samples agree 
within 2a. It is also particularly reassuring since Sgr 
members in the spectroscopic sample have very different 
colors (g— i ~ 0.75) from the Sgr members in the photo- 
metric sample (g— i ~ 0.3) of the Sgr giants. Therefore 
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the agreement of the proper motions from two samples is 
a proof of a small level of color-related systematic effects. 

The proper motion signals we have measured have a 
very large contribution of the solar reflex motion. Since 
we have distance estimate s to the Sgr stream in the South 
measured elsewhere (e.g. iKoposov et all 120121 ). we can 
correct for this and check whether the proper motions 
are actually aligned with the streams. Figure [5] shows 
the proper motion vectors after applying the solar reflex 
corrections. As we can clearly see, the proper motions 
are indeed properly aligned with the streams and are 
consistent with each other. For this calculation we have 
still used the distance to t he faint stream as given by 
IKoposov et al.1 (120121 12013I ). although, there is evidence 



m 
3- 



Slater et al.1 (|2013l) that the stream the faint stream is 



5 kpc closer. 



5. COMPARISONS TO EARLIER WORK 



There are many models of the Sgr stream in the lit- 
erature, all purporting to provide the distances, veloci- 
ties and proper motions as a function of position on the 
sky (e.g. iFellhauer et all 120061: iPenarrubia et all l2010t 
lLaw fc Maiewskil l2010al ). Although quite detailed, all 
the models fail to reproduce at l east some of the features 
that we see on the sky (e.g., iNiederste-Ostholt et al.1 
120101 : IKoposov et"al~ll2012t ). But, it still instructive to see 
where our data measurements lie relative to the existing 
models. We have chosen the lLaw fc Maiews ki (2010a) 
model as a comparison benchmark, because it is arguably 
the most comprehensive and up-to-date. 

There are several caveats to be borne in mind. The 
first is related to the fact that while our central mea- 
surement in Stripe 82 at a ~ 35° or (A, B) ~ (100°, 0°) 
corresponds directly to the center of the trailing tail in 
the simulation, the fainter stream which crosses Stripe 
82 at a ~ 15° or A,fl - (90°, 1 0°) doesn't have 
a coun terpart in the simulation by lLaw fc: Majewskl 
(|2010af ). The second is related to the choice of rota- 
ti on velocity of the L ocal Standard of Rest. As shown 
in iCarlin et al.l (|2012[ ). the observed proper motion sig- 
nal is sensitive t o the a dopted Vlsr- The models of 
lLaw fc Maiews ki (2010a) have been computed and fit- 
ted under the assumption of the IAU standard value of 
the V[ J sR=220kms~ 1 , while the curre nt best estimates 
are slightly higher at 230 - 250kms~ 1 (lB^yxei^|2009|; 
iReid et all 120091: iDeason et all 120111) . Throughout the 
paper, we adopted the value of 235 km s _1 . 

Fig. [6] shows our data po ints overplotted onto the 
distribution of tracers from lLaw &: Maiewskil (|2010al) 
model, where we have selected tracers from the trail- 
ing tail only (Lmflag = — 1 and Pcol <= 7), and lying 
within 20 degrees of Stripe 82 (\S\ < 20). The A values of 
our measurements correspond to the centers of the rect- 
angular areas used to perform the measurements, and 
the error-bars indicate the extent of those areas. The 
agreement for the central field (A ~ 105°), as well as 
for the edge of the bright stream (A ~ 120°), is quite 
good. At the location of the fainter Sgr stream (empty 
circle on the plot), the model doesn't have many par- 
ticles, as expected. In order to facilitate further com- 
parisons, Table [3] gives the positions and velocities cor- 
responding to our measurements. Here, (X, Y, Z) and 
(U, V, W) are in a right handed Galactocentric coordi- 
nate system with the Sun located at X = —8.5 kpc, and 
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Figure 6. The comparison of the measured Stripe 82 proper mo- 
tions (red w i th err or bars) with both the earlier measurements of 
Carlin et al. (2012, grey with error bars) and the simulations of 
lLaw fc Maiewski (2010a, blue, green, yellow dots). Our measure- 
ment of the proper motion of the fainter of the two trailing tails 
is id entified by an empty red circle. Note that the datapoints 
from Law & Maiewski (2010a) are selected from trailing tail only 
and from the region near Stripe 82 \S\ < 20°. The points are col- 
ored according to S such that points with 6 ~ —20° are dark blue, 
points with <5 ~ 20° are orang e, and points near S tripe 82 are light- 
green. Two of the fields from Carlin et al. (2012) lie within Stripe 
82, while two other fields are ~ 15 degrees away from the Stripe 
82. 

Vl sr = 235 km s" 1 . The distances taken from the work 
bv IKoposov et all (|2012f ). and the error on the velocities 
does not take into account any possible systematic error 
distance of ~ 10%. Readers wishing to compare Sgr 
disruption models with our observations should use Ta- 
ble [3] only for rough checks. A proper comparison entails 
measuring the average proper motions within the same 
area of the sky as we do and comparing these numbers 
directly with our Table [2j 

We also compa r e our measurements with the data 
from ICarlin et all (|2012f) . Out of 6 fields analysed by 
ICarlin et all (|2012| ). 4 have measurable Sgr signal and 
only 2 of those (SA 93 and SA 94) lie near the Stripe 82. 
The other 2 fields with detectable Sgr signal are ~ 15 
degrees away from the Stripe 82. As Figure [6] shows, the 
error-bars of our m easurements are much tighter than 
ICarlin et all ([2012), and the proper motions themselves 
agree approximately within the combined errors. 

6. CONCLUSIONS 

We have measured the proper motion of the Sgr stream 
using Stripe 82 data. By tying the astrometry to the 
known QSOs and combining the measurements for large 
samples of stars, we have been able to achieve high pre- 
cision proper motions. Our measurements have been 
performed for two distinct groups - spectroscopically se- 
lected red giants/subgiants with SDSS radial velocities 
and photometrically selected MSTO stars. The proper 
motions of those agree very well and after correcting for 
the solar reflex motion are tightly aligned with the Sgr 
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stream. Our results are in agree ment with earlier mea- 
surements by iCarlin et al.l (|2012l) . but, by virtue of the 
large numbers of stars in our samples, our statistical er- 
ror bars are substantially smaller, typically about 0.1 mas 
yr" 1 . 

There are three fields in Stripe 82 for which the proper 
motion of the photometric sample has been measured, 
and one field for the spectroscopic sa mple. Combin- 
ing th is information with distances from iKoposov et al.l 
(2012) gives us the full six dimensional phase space coor- 
dinates of the Sgr trailing stream at four locations along 
Stripe 82. We provide a table of three dimensional po- 
sitions and velocities of Sgr stream stars, in which the 
contributions of the motion of the Sun and LSR have 
been removed. 

To complement the work on the Sgr streams in the 
South, it would be particularly useful to carry out a cor- 
responding kinematical analysis for the streams in the 
North. This is a subject to which we plan to return 
in a later contribution. The combination of the proper 
motions, radial velocities and distances in both Galac- 
tic hemispheres should allow us to make further progress 
in the understanding of the complicated structure of the 
Sgr streams and solve some of the riddles posed by their 
existence. 
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